Theory of Cross Phase Modulation for the Vibrational Modes of Trapped Ions 
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We analyze nonlinear coupling between individual vibrational quanta for trapped ions. The 
nonlinear Coulomb interaction causes a Kerr-type Hamiltonian, for which we derive an analytical 
expression for the coupling constant x- In contrast to a previously published formula [T], our result 
is in close agreement with experimental data. 



It is well known that the nonhnear interaction at sin- 
gle quantum level is crucial to the development of quan- 
tum technologies. Systems in which two field modes are 
coupled by a Kerr nonlinearity (ie. with a Hamiltonian 
Hkerr = TT-xfirfia) enable quantum nondemolition mea- 
surements [2 , allow quantum gate operations in photonic 
quantum computation [3], [4], and assist nondestructive 
Bell-state detection [5]. However, it is often difficult to 
obtain the strong nonlinearities required to observe these 
effects at the single quantum level [6j,|7j. One exception 
is provided by the joint vibrational modes of trapped 
ions, which have an intrinsic nonlinearity considerably 
stronger than that of photons, and can be utilized to our 
advantage [5]. In this paper, we shall present an analyt- 
ical expression to describe this phonon-phonon interac- 
tion for single quanta and compare it to the experimental 
results. 

Let us concentrate on the dispersive cross-Kerr effects 
causing shifts of the normal mode frequencies [9] rather 
than the resonant mode-mode coupling [TU] . Suppose we 
have two equally charged ions of the same mass, confined 
in a three dimensional harmonic trapping potential. We 
assume the strength of the trapping potential is the same 
in two transverse directions, and is characterized by an 
angular frequency uj_\_, while in the axial direction, the 
trapping potential characterized by angular frequency lo^ 
is considerably weaker (Figure 1). We have the following 
classical expression for the system's kinetic energy and 
potential. 
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Figure 1: Schematic diagram of the system. 
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where the center of mass and the relative motions are 
written as X — 5(2:1 + X2), Y — ^{yi + 1/2), Z = 5(21 + 
Z2), and X = i(xi - X2), y = ^{yi - 2/2), z = ^{zi - Z2) 
respectively. Further, m is the mass of one of the ions, 
q is the charge of one of the ions (for singly ionized ions 
q = e, the positively valued fundamental charge), and eg 
is the permittivity of free space. 

Since the center of mass motion, characterized by 
{X, Y, Z} is entirely decoupled from the relative motional 
degrees of freedom {x, y, z}, we ignore it; this leaves us 
the following Lagrangian for the relative motion: 



'Electronic address: 'rebecca.nie@utoronto.ca' 
tEIectronic address: christian. roosSSuibk.ac. at 
J Electronic address: jdt'vj@pliysics.utoronto.ca| 



/• • 2 I -2 I ■2\ 

= m{x + y + z ) 



2 „2 



mujl^x +y ) + mu!^z' 



2 



where 



Defining the canonical momenta in the standard way (i.e. 
Px = dCrei/di, and similarly for py and Pz), this La- 
grangian gives us the Hamiltonian for the relative motion 
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Note that the effective mass for the relative motion is 
2m, hence the cancellation of the factor 2 in the de- 
nominator in the potential energy, and the appearance 
of 4 in the denominator of the kinetic energy. Rather 
than confuse things by introducing another symbol for 
this effective mass, we retain m, which is the mass of 
one of the the two ions. The equilibrium separation 
of the ions is determined in the usual way by solving 
dVrei/dx =dVrei/dy =dVrei/dz = 0. Provided lj± > oj^, 
the equilibrium separation of the ions is x = 0,y = 0, z = 
zo, where 
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Here a is the fine structure constant (and we have as- 
sumed the ions are singly ionized, i.e. q = e), ?i is 
Planck's constant, and c the speed of light. We shall 
write z = u + Zq, u being a new variable indicating the 
ions' displacement from their equilibrium separation. 

It is insightful to expand the potential Vrei around the 
equilibrium position (0,0, Zq)- First of all, let us write 
the potential in terms of an expansion constant ^: 

V[(0, 0, Zo) + <r(a;, y, u)] = Smw^Zg 
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The first term on the right hand side, which is indepen- 
dent of represents the constant potential energy of the 
equilibrium, and has no effect on the dynamics. The 
terms 0(<;^) are harmonic potentials: immediately, we 
see from this term that the natural angular frequency of 
the "ro cking mo de" oscillations in the x, y-directions is 
Ur = — i^z, while the "stretch mode" oscillation in 

the z direction has the angular frequency LUg = VSuj^- 
Keeping terms of 0('?'*) and setting the expansion con- 
stant c = 1, the Hamiltonian becomes 
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represents the harmonic motion (which, to a very good 
approximation, is the dominant feature of the ions' mo- 
tion); 
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and 



are two terms which represent the lowest order perturba- 
tions of the harmonic terms. 

To quantize the motion, we introduce the following 
operators: 
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The quantized Hamiltonian is 

H = Ho + Hi, 

where 



{c + c^), Pz = i\/m?i(Js{c - c^). 



and 



Ho = hu!r{nx + 'hy + l) + hws{hs + -j), 



Hj = V(^) + V(^). 



fix, fly, and hs are the phonon number operators for their 
respective normal modes. 

Now we are in a position to treat Hi with the standard 
second-order perturbation theory, but before we begin, 
recall that our goal is to extract the x term from the Kerr- 
type Hamiltonian Hkerr = Tixfirfis- To do that, we would 
like to calculate the shift in the stretch mode frequency 
as an effect of the coupling, when one phonon is inserted 
into the stretch mode. In other words, we want 

^ e(n^ + 1, n%, nf) - e{ns,rv^, nv) ^^^^ 
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Let us start with t^*^^^ , given by eq.([9|; which, when quan- 
tized, becomes 
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Here ^ = {2fujL!z/ oi^rac^Y^^ is a dimensionless term of the 
order of 10~^ for atomic ions. The off resonance terms, as 



well as the second term in Eq.( 11 1 will not have any effect 
to first order on the energy. Further, we shall neglect the 
third term, as it involves anharmonicity of the stretch 
mode rather than cross-coupling between the stretch and 
rocking motion. Thus we find that V^'^'^ has the following 
contribution to the energy shift: 
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Now we deal with V'^'^'^ , eq. m\ , which can be quantized 
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where C = y C^^'^f^z/32ci;^ is a constant with units of 
energy. 

Some facts about the perturbation V'^^^ are now ap- 
parent. Firstly, V'^^^ represents the creation and anni- 
hilation of vibrational quanta, and does not yield any 
value change in energy to first order. Therefore, to cal- 
culate energy shifts, we need to seek out the cross-coupled 
terms from in its second order perturbation. Secondly, 
since the terms involving (? and c^^ do not give rise to 
any cross-coupled terms, they can be safely neglected. 
Finally, only a few states out of the infinite number of 
possible states can yield nonzero value during the second 
order perturbation of the rest of V'^'^\ All those states, 
their corresponding matrix elements, and their respec- 
tive differences between the unperturbed Hamiltonian's 
eigenvalues are listed in Table 1. 
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Table 1: The states \<fim) = |m^,m^,ms) relevant to the 
second order perturbation of V^^\ their corresponding matrix 
elements {{(fin\V'-^''\(fim} = (n^ , n^, n^l F'^^lm^, mjf, ttis)), and 
their respective differences between the unperturbed Hamilto- 
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parameter = ^y (Ti'^u)'iu}z/'i2ujr is a constant with units of 
energy. 



The next step is to sum these matrix elements in the 
standard manner given by second order perturbation the- 
ory , i.e.. 



A£;: 



E l (n^,ny, ns\V^^^ \mf, m^, Ws)|^ 



where {m^,m^,ms}' denotes the condition m^myrris ^ 
n^n'^ns [H]. Keeping only the cross-coupled terms, we 
find that the energy shift given by F'-'^-' is 

e(3)«,n5f,n,) = 

{n^,+nv + l). (14) 
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Combining Eq.( 12 1 and Eq.( 14 1, the energy shift is found 
to be 
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Higher order perturbations, obtained by keeping terms 
of higher order in <; in eq. ([s]) , give rise to terms in higher 
powers of ^, which, as mentioned above, is a small quan- 
tity for atomic ions. Thus the contributions of these 



higher terms may safely be neglected. From Eq.(15), we 
find the expression for the cross-Kerr coefficient x 
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This formula differs from eq.(l) of ref.[T] in that the first 
term in the square bracket is 1/2 rather than 1. The 
distinction arises entirely from the more thorough treat- 
ment of the second-order shift in energies due to the V-^^ 
presented in this paper. 
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Figure 2: Shift of the stretch mode frequency by a single 
phonon as a function of the axial trap frequency u^. The 
solid line is the theoretical model predicted by Eq.(16l, and 
the points are the data. 



Our theoretical model can be readily compared to 
the experimental data presented by Roos et ai., where 
two ''^Ca^ ions are confined in a linear Paul trap [T]. 
Their experimental results for the measured cross-Kerr 
coupling constant as a function of lOz are compared to 
Eq.(16), and presented in Figure 2. As we can see, the 



theoretical model and the experimental results are in 
good agreement, much better than the theoretical curve 
presented in ref. [I:- 

In conclusion, we have derived an analytical expres- 
sion for the cross-Kerr coupling constant of two equally 
charged ions of the same mass trapped in a an-axial har- 
monic potential, while showing that the theoretical model 
is verifiable via experimental means. Since the analytical 
expression is reasonably general, it may describe other 
trapped ion systems giving rise to a cross coupling of 
harmonic oscillators governed by a Kerr-like Hamiltonian 
H cx fisfir- On the other hand, it would be interesting to 
further investigate the reason behind the experimental 
deviation from the theoretical model in high u)z regime. 
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